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Foreword 

The  following  is  our  final  report  on  the  project  entitled  "Electromagnetic  model¬ 
ing  of  quantum- well  infrared  photodetectors,"  which  started  on  20  April,  1999,  and  was 
completed  (after  being  extended)  on  August  19,  2003. 

The  principal  aim  of  this  research  work  was  to  develop  analytical  models  of  peri¬ 
odic  configurations  that  occur  mostly  in  quantum-well  infrared  photodetectors  (QWIPs), 
as  well  as  in  other  optoelectronic  devices.  To  strengthen  this  theoretical  study,  a  close 
collaboration  was  maintained  with  the  Dr.  K.-K.  Choi’s  group  at  the  Electro-Optics  and 
Photonics  Division,  Army  Research  Laboratory,  Adelphi,  MD,  and  with  Prof.  D.  C. 
Tsui's  group  at  the  Department  of  Electrical  Engineering,  Princeton  University,  Prince¬ 
ton,  NJ.  This  has  enabled  us  to  apply  our  electromagnetic  modeling  techniques  to  a  vari¬ 
ety  of  QWIP  configurations  under  active  experimental  investigation  by  those  two  groups. 
We  have  thus  arrived  at  interesting  and  useful  explanations  of  their  measurements  and 
have  provided  guidance  in  designing  their  experimental  efforts. 

The  major  achievement  of  our  study  was  the  development  of  a  rigorous  modal 
approach  for  the  analysis  and  design  of  multilayered  QWIP  configurations  that  incorpo¬ 
rate  arbitrarily  shaped  periodic  contours,  and  contain  materials  consisting  of  lossy  dielec¬ 
tric  media  and  metallic  conductors.  This  modal  approach  provides  physical  insight  into 
the  detection  process,  on  the  one  hand,  and  readily  perfonns  an  accurate  evaluation  of  the 
photocurrent  characteristics,  on  the  other  hand.  Furthermore,  our  approach  employs 
transmission-line  techniques  that  serve  as  a  very  powerful  tool  in  the  engineering  design 
of  QWIP  structures  having  specialized  and/or  desirable  characteristics,  such  as  high  effi¬ 
ciency,  two-color  detection  capability,  narrow  or  wide  bandwidths,  etc.  As  described 
herein  the  theoretical  basis  of  our  rigorous  analytical  approach,  its  associated  transmis¬ 
sion-line  and  computational  techniques,  and  its  application  to  a  wide  variety  of  QWIP 
configurations  have  been  published  in  the  professional  literature  and  reported  at  both 
general  and  specialized  professional  meetings. 


Statement  of  the  problem  studied 


Metal 


Substrate 


Quantum  well  infrared  photodetectors  (QWIPs)  generally  consist  of  a  quantum 
well  (QW)  region  inserted  between  conducting  layers  that  serve  as  electrodes.  Additional 
layers  are  also  present  due  to  specific  device  requirements  or  fabrication  purposes.  The 
material  forming  a  QWIP  are  GaAs  and  AlxGai_TAs,  but  other  materials  (such  as  MgF2, 

epoxy,  gold  or  other  conductors) 


are  also  part  of  the  complete 
QWIP  device.  Because  only  the 
incident  electric  component  nor¬ 
mal  to  the  QW  layers  generates 
photocurrent,  it  is  necessary  to 
incorporate  a  mechanism  that  de¬ 
flects  the  field  into  that  desired 
direction.  This  can  be  achieved 
by  a  grating  configuration,  which 
can  be  above  or  below  the  QW 
layer.1"3  Alternatively,  the  QW 
layer  can  itself  be  fonned  into  a 
periodic  structure,  thus  serving  as 
the  required  grating.4,  5  Because 
the  latter  situation  is  easier  to  fab- 


Trapezoidal 


Rectangular 


Undercut 


Fig.  1:  Typical  lamellar  QWIPs 


ricate  and  exhibits  a  variety  of  other  advantages,  we  shall  assume  that  this  is  the  case  for 
illustration  purposes. 

Several  QWIP  configurations  are  shown  in  Fig.  1,  where  the  periodic  portions 
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Air  or  epoxy  cover  i 


Air  or  GaAs  substrate 

(b)  yz  cross-section 

Fig.  2:  Grid  QWIP  with  gold  and  MgF2  layers 

consist  of  rectangular  or  trapezoidal  cells,  which  may  be  symmetric  or  slanted,  either  in- 
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tentionally  or  due  to  the  fabrication  (etching)  process.  Because  of  their  basic  lamellar 
form,  these  QWIPs  can  be  treated  by  a  simple  two-dimensional  analysis.  A  more  com¬ 
plex  QWIP  structure  is  shown  in  Fig.  2,  where  all  dimensions  are  in  pm.  A  large  number 
of  layers  and  different  materials  are  present  and  the  materials  include  gold  (Au),  whose 
electromagnetic  behavior  is  very  much  different  from  that  of  the  dielectric  materials  oc¬ 
cupying  most  of  the  QWIP  configuration.  Because  the  QWIP  in  Fig.  2  is  periodic  along 
both  horizontal  directions,  it  fonns  a  grid  having  trapezoidal  xz  and  hourglass  yz  cross- 
sections,  which  involves  a  three-dimensional  geometry.  Obviously,  finding  the  electro¬ 
magnetic  fields  due  to  energy  incident  on  such  a  configuration  poses  a  formidable  bound¬ 
ary-value  problem  if  an  accurate  electromagnetic  field  solution  is  required. 

To  accommodate  the  large 
variety  of  structures  illustrated  in 
Figs.  1  and  2,  we  have  explored6' 7  the 
general  situation  described  in  Fig.  3, 
where  an  electromagnetic  field  is  in¬ 
cident  from  below  into  the  substrate 
of  a  multilayered  stack  containing 
periodic  elements  having  arbitrary 
shapes  and  materials.  The  configura¬ 
tion  in  Fig.  3  allows  for  arbitrarily 
shaped  periodic  regions,  and  the  ma¬ 
terials  can  generally  be  lossy  biaxial 
dielectrics  and/or  metallic  media,  in¬ 
cluding  perfect  electric  conductors 
(PECs).  As  shown  here,  the  geome¬ 
tries  in  Figs.  1  and  3  are  two- 
dimensional  (2D)  and  refer  to  lamel¬ 
lar  QWIP  configurations.  However, 
E  we  have  also  extended  our  results  to 

Fig.  3:  General  geomehy  of  QWIP  configurations  three-dimensional  (3D)  grid  configu¬ 
rations,  such  as  that  shown  in  Fig.  2, 
by  allowing  for  periodic  variations  along  the  y  coordinate  in  Fig.  3. 

In  all  cases,  we  have  developed  rigorous  analytical  solutions  that,  as  described  in 
the  next  Section,  yield  the  electromagnetic  fields  at  any  point  in  the  QWIP  structure.  By 
obtaining  these  fields,  it  is  then  possible  to  derive  the  power  absorption  in  the  QW  layers 
and  thus  evaluate  the  photocurrent.  This,  in  turn,  detennines  the  quantum  efficiency. 
Other  characteristics  of  the  C-QWIP  device  can  similarly  be  obtained  from  the  field 
quantities.  We  have  applied  these  analytical  tools  to  construct  computational  programs 
and  obtained  results  for  the  experimental  QWIPs  developed  by  the  ARL  and  Princeton 
University  groups.  As  described  further  below,  the  agreement  between  our  theoretical 
modeling  and  their  experimental  data  was  very  high. 
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Summary  of  the  most  important  results 

The  results  achieved  by  our  work  fall  into  two  principal  categories.  The  first  one 
is  the  development  of  a  rigorous  analytical  solution  for  the  electromagnetic  problem 
posed  by  the  general  situation  shown  in  Fig.  3.  The  second  category  is  the  application  of 
this  solution  to  a  wide  variety  of  actual  QWIP  structures.  Of  these,  most  have  already 
been  (or  are  expected  to  be)  fabricated  and  tested  by  the  ARL  and  Princeton  University 
groups.  In  this  context,  we  have  developed  criteria  for  designing  QWIP  configurations 
having  novel  functions  or  desirable  improved  characteristics.  Based  on  these  design  con¬ 
siderations,  new  QWIP  configurations  are  either  being  fabricated  for  immediate  testing, 
or  under  consideration  for  future  applications,  as  further  described  below. 

A.  Rigorous  electromagnetic  modeling 

The  requisite  analytical  models  for  the  general  periodic  structure  in  Fig.  3  were 
developed  sequentially  as  follows: 

(i)  It  was  first  assumed  that  the  materials  were  dielectrics  having  complex  refractive 
index,  which  included  metals  but  not  perfect  electric  conductors  (PECs).  Periodic  re¬ 
gions  having  arbitrary  profiles  were  subdivided  into  sufficiently  many  j  =  1,2,  ...  J 
thin  layers,  each  of  which  was  then  approximated  by  a  step-wise  periodic  region  of 
height  tj  and  uniform  periodic  regions  having  dielectric  constants  Sj,  sr\,  cri,  ...  ,  as 

shown  in  the  periodic  layer  of  Fig. 
3.  By  using  a  rigorous  modal  ap¬ 
proach,  the  scattering  of  a  plane 
wave  by  such  a  structure  can  then 
be  described  in  terms  of  an  equi¬ 
valent  transmission-line  network.6 
Specifically,  for  a  relatively  simple 
situation  having  a  periodic  layer 
embedded  in  regions  having  ho¬ 
mogeneous  dielectric  media,  such 
a  network  is  shown  in  Fig.  4. 

The  transmission  lines  shown 
here  are  given  by  matrix  forms  of 
characteristic  admittances  Y )  and 
propagation  factors  K 7,  whose  in¬ 
dividual  terms  represent  the  vari¬ 
ous  diffracted  orders  supported  by 
the  periodic  region.  All  of  these 
lines  are  connected  in  tandem,  except  that  interface  transfonners  are  needed  at  the 
junctions  with  the  periodic  layer.  By  using  units  such  as  that  shown  in  Fig.  4,  it  is 
possible  to  construct  the  solution  to  any  arbitrary  geometry  as  in  Fig.  3. 

At  this  stage  of  our  project,  the  analysis  was  restricted  to  2D  geometries  of  the 
lamellar  form,  but  the  dielectric  media  could  be  uniaxial,  so  that  QW  regions  were 
easily  accommodated.  For  these  cases,  the  interface  transformers  shown  in  Fig.  4  had 
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Fig.  4  -  Typical  portion  of  a  grating  geometry 
and  its  equivalent  network  representation 
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simple  matrix  relations,  which  were  determined  by  requiring  that  the  fields  of  each 
diffracted  order  were  continuous  at  the  interface. 

(ii)  The  above  phase  provided  a  convenient  modeling  capability  for  devices  contain¬ 
ing  dielectric  media,  as  well  as  metals  in  the  form  of  thin  horizontal  plates.  However, 
for  components  consisting  of  highly  conducting  metals  (e.g.,  gold  or  silver)  having 
thicker  dimensions,  the  corresponding  computational  program  converged  slowly.  To 
avoid  this  aspect,  the  second  phase  of  our  study7  extended  the  analysis  to  geometries 
that  included  PECs,  which  were  found  to  provide  excellent  numerical  results  for 

situations  involving  such  highly  conducting  media.  The 
presence  of  such  PECs  does  not  change  the  transmission¬ 
line  aspect  depicted  in  Fig.  4,  but  it  strongly  affects  the 
analytical  expressions  for  the  interface  transformers. 

Some  of  the  geometries  involving  PEC  materials  are 
shown  in  Fig.  5,  for  which  the  interface  transformations 
have  been  analytically  derived  and  applied  to  specific  prac¬ 
tical  situations.  In  conjunction  with  the  transmission-line 
networks  illustrated  in  Fig.  4,  these  transformations  served 
as  templates  for  computational  programs  that  provided  ac¬ 
curate  results  for  the  operation  of  actual  lamellar  devices 
whose  perfonnance  had  been  experimentally  determined. 
Because  the  agreement  between  the  experimental  data  and 
the  theoretical  model  was  very  good  to  excellent,9’  10  our 
modal  approach  can  readily  serve  as  a  powerful  tool  for 
both  analysis  and  design  purposes.11'13 

(iii)  The  preceding  discussion  has  dealt  primarily  with  2D 
lamellar  configurations.  For  the  more  general  situation 
having  periodic  variations  with  respect  to  both  x  and  y  co¬ 
ordinates,  e.g.,  for  QWIPs  of  the  grid  type  in  a  focal  plane 
array,  a  3D  geometry  is  needed.  We  have  therefore  ex- 
tended  our  2D  modal  transmission-line  approach  to  3D 
and  have  successfully  examined  several  practical  situa¬ 
tions.  In  particular,  the  complex  geometry  shown  in  Fig.  2  was  examined  and  design 
criteria  were  established  for  particular  focal-plane  arrays.  Such  arrays  have  been 
fabricated  and  preliminary  data  have  shown  good  agreement  with  our  theoretical  re¬ 
sults,  but  the  work  is  still  ongoing  and  a  more  thorough  comparison  between  our 
modeling  approach  and  the  actual  experimental  data  remains  to  be  obtained  at  a  later 
date. 


(a)  Metal-dielectric  interface 


(c)  Offset  metal-metal  interface 
Fig.  5:  Typical  metallic  interfaces 


B.  Evaluation  and  design  of  specific  QWIP  configurations 

The  rigorous  modeling  approach  described  above  was  applied  to  evaluate  the  per¬ 
formance  of  experimental  QWIPs  and  to  design  new  configurations,  as  described  below. 

(i)  The  performance  of  lamellar  QWIPs  having  a  slanted  trapezoidal  profile  (as  shown 
in  Fig.  1)  was  evaluated  and  the  theoretical  results6  showed  excellent  agreement  with 
the  experimental  data.9  This  agreement  was  especially  significant  since  the  QWIP 
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performance  was  analytically  predicted  without  any  adjustable  parameters,  on  the  one 
hand,  and  its  geometry  involved  an  asymmetric  grating,  on  the  other  hand. 

(ii)  Design  criteria  were  established  for  lamellar  QWIPs  having  rectangular  gratings 
in  a  grid  configuration.10  The  top  of  the  rectangular  corrugations  were  covered  by  a 
thin  gold  layer.  As  in  item  (i)  above,  the  agreement  between  the  analytical  predic¬ 
tions  and  the  experimental  results  was  excellent.  Most  importantly,  we  have  thus 
shown  that  the  metal  top  dimensions  play  a  critical  role  in  maximizing  the  photocur¬ 
rent.  Specifically,  that  current  peaks  if  the  width  of  the  metal  layer  is  equal  to  an  odd 
number  of  half  wavelengths  in  the  material.  This  occurs  because  the  induced  current 
in  the  metal  acts  as  an  electric  dipole  whose  behavior  serves  as  an  effective  design 
criterion  to  optimize  performance  at  desirable  wavelengths. 

(iii)  Based  on  the  results  obtained  in  item  (ii)  above,  we  have  established  the  ground- 
work  for  constructing  a  spectrometer  using  QWIP  elements.  ’  '  This  device  consists 
of  channels  containing  those  elements  in  a  grid  pattern.  The  dimensions  of  the 
QWIPs  in  the  channels  are  tailored  to  provide  response  peaks  at  different  wave¬ 
lengths.  Design  criteria  were  developed  for  a  spectrometer  containing  up  to  10  such 
channels  and  operating  in  the  range  7.5  -  12.0  pm,  as  well  as  in  a  broader  6.0  -  15.0 
pm  range.  Preliminary  experimental  results  support  this  concept  and  its  design  feasi¬ 
bility. 

(iv)  A  broad  variety  of  QWIPs  having  triangular  or  trapezoidal  grating  profiles  were 
examined  '  ’  under  different  conditions,  e.g.,  with  or  without  metal  tops,  with  air 
or  epoxy  covers,  placed  over  thick  or  thinned  substrates,  etc.  The  results  have  shown 
that  the  placement  and  dimensions  of  the  metal  tops  may  enhance  performance  in  cer¬ 
tain  cases,  but  the  presence  of  epoxy  cover  and/or  additional  insulating  layers  (usu¬ 
ally,  MgF2  material)  could  degrade  performance.  All  of  these  considerations  have 
provided  a  plethora  of  design  considerations,  most  of  which  have  been  verified  ex¬ 
perimentally. 

(v)  The  complex  grid  structure  shown  in  Fig.  2  for  constructing  a  1024  x  1024  focal 
plane  array  was  analyzed  and  numerical  results  have  shown  that  it  can  provide  an  ef- 
fective  performance  over  a  broad  7.5  -  11.5  pm  wavelength  range.  ~  This  desirable 
performance  is  satisfactorily  maintained  if  the  structure  uses  an  epoxy  cover  and  has  a 
thin  substrate.  Hence  such  a  focal  plane  array  can  serve  as  the  active  component  in 
an  infrared  imaging  camera. 

(vi)  Using  3D  analysis,  we  have  evaluated  the  performance  of  grid  QWIPs  in  which 
each  lattice  contains  a  post  incorporating  QW  layers.  The  posts  may  have  circular, 
square  or  triangular  cross-sections,  thus  leading  to  different  detection  efficiencies, 
which  are  subject  to  the  symmetry  posed  by  the  relevant  cross-sections.  The  results 
show  that  the  circular  posts  provide  optimum  performance  and  their  high  symmetry 
renders  them  least  susceptible  to  variation  in  fabrication. 
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